MCTP-11-37 



D vs d: CP Violation in Beta Decay 
and Electric Dipole Moments 

John Ng 

TRIUMF Theory Group, 4OO4 Wesbrook Mall, 
Vancouver, BC, Canada, V6T 2A3 

Sean Tulin 

Michigan Center for Theoretical Physics, Department of Physics, 
University of Michigan, Ann Arbor, MI, USA, 48109 
(Dated: November 4, 2011) 

Abstract 

The T-odd correlation coefficient D in nuclear /3-decay probes CP violation in many theories 
beyond the Standard Model. We provide an analysis for how large D can be in light of constraints 
from electric dipole moment (EDM) searches. We argue that the neutron EDM dn currently 
provides the strongest constraint on D, which is 10 — 10^ times stronger than current direct limits 
on D (depending on the model). In particular, contributions to D in leptoquark models (previously 
regarded as "EDM safe") are more constrained than previously thought. Bounds on D can be 
weakened only by fine-tuned cancellations or if theoretical uncertainties are larger than estimated 
in dn- We also study implications for D from mercury and deuteron EDMs. 
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I. INTRODUCTION 



The search for CP violation beyond the Standard Model (SM) remains an open question 
at the forefront of nuclear physics, particle physics, and cosmology.^ New CP violation is 
a generic feature of physics beyond the SM and is likely required to explain the baryon 
asymmetry of the Universe j2|. Furthermore, unlike the SM Kobayashi Maskawa (KM) 
phase [ij, new CP violation may be unconnected with flavor and can be probed in systems 
of "ordinary matter" through searches for T violation in nuclear /3-decay and electric dipole 
moments (EDMs) of atoms, nucleons, and nuclei. 

CP violation in /3-decay is manifested through T-odd triple product correlations Q • (See 
Refs. jH, IgI 0, [sj for reviews of fundamental symmetry tests in /3-decay.) In this work 



we 



study the so-called D correlation, corresponding to the triple product (J) ■ pgXp,^, where 
(J) is nuclear polarization, and Pe (pi^) is the {v) momentum. It is useful to write 
D = Df + Df to delineate fundamental T violation (Dt) from T-even final state effects 
(Df) In the SM, the KM phase contribution to Df is vanishingly small jof. Therefore, 
to the extent that Df is computable or negligible, measurements of D directly probe CP 
violation beyond the SM. 

To date, D has been measured for the neutron hiA hjj, 12, 13, 14, uM and ^^Ne 



16U171 



The best neutron D measurement has been obtained recently by the emiT collaboration [15 



-1.0 ±2.1) X 10' 



(1) 



Final state interactions give Df = O{10~^) [l8|, and have been computed to an accuracy 
better than 1% [l^. Although Dn measurements so far agree with SM expectations, there 
remains (in principle) a discovery window for future experiments down to D„ ~ 10~^. For 
^^Ne, an average of previous measurements 16|, llTf gives 



D 



Ne 



;i ±6) X 10- 



(2) 



which has reached a level comparable to final state interaction effects ~ 10 [Ti 

Measurements of EDMs (denoted d) are also sensitive to CP violation in and beyond the 



SM (20| . No EDM has yet been observed, but many future experiments await [2l| . Currentb 
the most significant EDM bounds are for the neutron 22 1, atomic mercury (^^^Hg) |23| 
atomic thallium (^"^Tl) jl^], and recently molecular YbF 



25 



atomic mercury 

These null results provide 
important constraints on CP violation in the SM due to the ^qcd phase associated with 
the strong interaction (present limits on dn require ^qcd < lO"^'^ |26|), and on CP violation 



beyond the SM, such as in the Minimal Supersymmetric Standard Model (MSSM) |27l . |28 
On the other hand, these observables are rather insensitive to the KM phase, requiring many 
orders of magnitude increases in sensitivities (see Ref. j20| and references therein). 

In this work, we compare D vs EDMs (in particular, dn and c/ng) as probes of CP 
violation beyond the SM. For a given model, any CP-odd phase contributing to D generates 
an "irreducible" EDM that can only be avoided by fine-tuned cancellations with other phases 
in the model. We compute the resulting bounds on D from EDMs in several new physics 
models: left-right symmetric models |29|, MSSM with R-parity violation 30|, models with 



^ The discrete symmetries discussed herein are charge conjugation (C), parity (P), and time reversal (T) 



symmetries. Assuming CPT invariance, T violation is equated with CP violat: 
^ We have added in quadrature statistical and systematic errors quoted in Ref. 



HQ. 
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exotic fermions [31i, and leptoquark (LQ) models [3^. Most of these scenarios, and the 
resulting constraints from EDMs, have been studied previously jsl, Iss], 34 1. Here, we provide 
several improvements: 



We take into account recent improved computations of dn [35| and dug |36| . Large un- 
certainties in the sensitivity of dug to the CP-odd isovector pion-nucleon coupling 36 
have weakened this constraint, and the dn bound currently provides the strongest limit 
on A- 

In the literature, LQ contributions to Dt are regarded as being safe from EDM con- 
straints (gI, 13] • We argue that Dt is in fact more constrained than previously thought. 



We also study implications for D from LQ searches at hadron colliders. 

• We compute for the first time Dt in the R-parity violating MSSM (with baryon-number 
violation), arising at one- loop order. 

• We provide a (partially) model-independent analysis that applies to all the aforemen- 
tioned models except LQs, for which the current limit on dn implies Dt < 3 x 10~^. 

We emphasize that D is much cleaner theoretically than the EDMs constraining it, which 
rely on hadronic and nuclear computations. Moreover, any realistic model may contain 
many different CP-odd phases, to which Dt and EDMs are sensitive to different linear 
combinations. The bounds we derive may be negated if there exist accidental cancellations 
between phases entering EDMs, and we neglect this possibility in our analysis. 

Our work is organized as follows. In Sec. [TTl we review CP violation in /3-decay. We 
also summarize theoretical computations of neutron, mercury, and deuteron EDMs from 
underlying CP- violating operators most relevant for constraining Dt. In Sees. IIIII and IIV[ 
we study constraints on Dt from EDM bounds in several scenarios beyond the SM, focusing 
in particular on LQ models. We present our conclusions in Sec. |Vl 



II. CP- VIOLATING OBSERVABLES 



A. Beta decay 

The most general set of /3-decay interactions can be parametrized at the quark level by 
an effective Lagrangian jsf 

= -^^^1^ ^ eaV^Pe uV^dp + h.c. (3) 

The chiralities (L, R) of the electron and down quark are labeled by a, /3. The index 7 = 
S", V^,T labels whether the interaction is scalar (P*^ = 1), vector (P^ = 7^^), or tensor 
(P^ = a^'^ /\/2). CP invariance is preserved in /3-decay if all ten complex coefficients 

'^LL5 '^L_R5 '^flL5 O.RR^ '^LL' '^LH' "^iJL' "^flij' '^LiJ' '^Sl (4) 

have a common phase (a|^/^, terms are identically zero). At leading order in the SM, all 
parameters vanish except a}[j^ = l. SM radiative corrections and new physics contributions 
to a\j^ can play an important role in the extraction of Vud (see, e.g., Refs. [13]), but for 
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CP-violating observables they can be neglected as subleading effects. We also hereafter set 
Vud = 1; correlations between D and EDMs depend on \Vud\, but the 0{few%) deviation 
from \ Vud\ = 1 is irrelevant compared to other theoretical uncertainties. We neglect possible 
flavor constraints by considering only couplings between first generation fermions. Lastly, 
we assume that /3-decay processes involve a single neutrino flavor eigenstate u^, and we 
allow for both L, R chiralities. Coefficients involving (sterile) right-handed neutrinos are 
only relevant provided these states are kinematically allowed _in /3-decay.^ 
In terms of the parametrization in Eq. ([3]) , D is given by 



Dt = Klm [alkali + a^^a^^) + k Im {al^al*^ + a|+a]J^) (5) 

gvQA 

where af_|_ = (af -|- af and aj^j^ = {a^^ + a^j^). For initial (final) state nucleus of spin J 
(J'), the coefficient k is 



^ AgygAMpMcT / J j 0.87 for n , . 

" ~ glMl + glMl^ M J + 1 ~ \ -1.03 for ^^Ne ' 



where gy = 1, (^a ~ 1-27 [40|, and Mp (Mqt) is the Fermi (Gamow- Teller) matrix element. 
Scalar and tensor form factors qs.r, originally estimated in Ref. j4l|, have been computed 



using lattice techniques (see Ref. j42|] and references therein). In this work, we neglect 
the scalar-tensor term in Eq. ([5]). The R coefficient, corresponding to the T-odd /3-decay 
correlation (J) ■ cTg x pe where CTg is polarization, has greater sensitivity to scalar- and 
tensor-type CP violation [sl, [6|. Moreover, these couplings are correlated with CP-odd 
tensor and scalar electron-nucleon coupliiigs, which are strongly constrained by ^^^Hg 2^ 
and 205T1 Q EDM bounds, respectively 0, H, Esl, H, . 



B. Electric Dipole Moments 

EDM searches are sensitive to a wide class of CP- violating operators that can arise beyond 
the SM: CP-odd quark and lepton dipole moments, Weinberg's three-gluon operator jisl ]. 
and four-fermion operators. Here, the most relevant one is a CP-odd four-quark operator 
Olr, given by 

4:Gf _ - - _ 

■^cff = — kiR Olr , Olr = i{uLl^dL dR-ff,UR - (IlI^ul UR-f^dR) (7) 

where kLR is the operator coefficient (normalized to 4:Gp/\/2). Within the context of left- 
right symmetric models, this effective interaction arises from CP-violating W-W mixing 
and has been studied previously jiol 50, 51, [Hi. We show in Fig. [1] that, by connecting 



the leptonic legs in a one-loop diagram, the same interference terms o^rO^I ^^l'^^r 
contributing to Dt also generate Olr- Moreover, this diagram does not involve any chirality- 
changing mass insertions, and therefore is not suppressed by any light fermion masses. Other 
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Sterile neutrinos with eV-scale mass have been studied reeently in eonncction with various neutrino 
anomalies (see, e.g., Refs. 38|), and important constraints are provided by cosmology 39|. We do not 



attempt to accommodate these issues here. 
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FIG. 1: CP violation entering Dt = 
operator Olr = i{uL7^dL (IrJ^iUr 
deuteron EDMs. 



K/m(a£^^a]^£ + a]^^a)j|j) automatically generates the four- quark 
- dL^^ULUR^^dpt), which contributes to neutron, mercury, and 



CP-odd operators (e.g., quark EDMs) also arise from new physics entering Z^j, but are 
suppressed by light masses and will not be considered here. 



The most significant EDM constraints on Olr are for the neutron |22[ and mercury 
atom 1231: 



< 2.9 X 10-^" ecm (90% CL) , I^Hgl < 3.1 x 10"^'' e cm (95% CL) 



(8) 



Future measurements of the deuteron EDM c//), expected at the level of 10 e cm or 



better (53| . will also provide important constraints on Olr- 

Ref. 35[ has performed a systematic computation of dn from CP-odd four-fermion oper- 
ators, using a combination of chiral effective theory techniques and quark model estimates 
for the hadronic matrix elements. Using their results, we take 



-1 X 10 



-19 



:lr e cm 



(9) 



with an 0{1) uncertainty on the numerical prefactor [35|. Earlier results |44j . l49l . l5lL [5 



are consistent at the order-of-magnitude level, but according to Ref. (35| are not as reliable 
in that they take into account different subsets of the full set of contributions to dn- 

Diamagnetic atoms (e.g., ^^^Hg) are also sensitive to CP-odd four-quark interactions. 



Interpretation of these measurements is a three step process (see, e.g., Ref. [4^, |45|). First, 
atomic calculations relate the measured EDM to the nuclear Schiff moment S- For the case 



of mercury, we take [55 



-2.6 X 10 ^''ecm X 



efm"^ 



(10) 



This numerical value (2.6) agrees with an earlier result (2.8) by two of those authors [56 



while another recent computation found a larger value (5.1) [57|. Second, the Schiff moment 



^ This value is consistent with a naive estimate dn ~ cMqcd / , where Mqcd ^ 1 GeV is the QCD scale 
and A is the scale of CP violation. Taking ~ GpkLR, we have ^ 2|fcLfi'| x 10~^^ ecm. Also, it is 
useful to note Olr = {uudi^^d— ui"fQudd + 6ut°'udij5t°'d— Qui^^t°'udt°'d) /i using a Ficrz transformation, 
where t"' is the 5'J7(3)c generator, to make contact with the notation of Ref. 35|. 
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is computed in terms of P, T-odd nucleon-pion couplings, of which only the isovector coupling 
gi is relevant since Olr is isovector pSl. Previous nuclear computations found (keeping only 
Qi terms): Sug = — 0.071 ggi ecm^ and Sug = — 0.055 ggi ecm^ ^o[, where g ~ 13.5 is 
the (CP-even) pion-nucleon strong coupling. However, a recent and improved computation 
by Ref. j36| found that the cji coefficient is very sensitive to the mo del- dependent nuclear 
potential inputs and may be suppressed by an order of magnitude (or more) and may have 

1591, H 



opposite sign compared to Refs. 
constraining Dt using dug 
remaining agnostic as to the sign (see Ref. 
Ref. jiol, we conservatively take gi = 2 x 



These nuclear physics uncertainties are crucial for 
In light of this unresolved issue, we take |5'Hg| = 0.0l5f|^i| efm^, 
for additional discussion). Third, following 
kLf>. Ref. j62| found a larger numerical 
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10" 

prefactor by a factor of 7. Putting all these pieces together, we take 



\d 



Hgl 



7x 10-24|fc 



lr\ ecm 



(11) 



with an uncertainty at the order-of-magnitude level. 

The deuteron EDM provides a much cleaner probe of gi compared to d yj^. Following the 
recent computation of Ref. j63| (in good agreement with earlier results jo^, 65, [g^), we take 

4.5 X 10"^° Ikr 



\d 



D\ 



1.9 X 10 



-14 I- I 

Ifl'il ecm 



iLR\ ecm 



(12) 



with (9(20 — 30%) uncertainty on the numerical factor (1.9) |63l. |67 



III. MODEL-INDEPENDENT BOUNDS ON D 

New physics contributions to /3-decay can be organized in terms of a hierarchy of non- 

— 1/2 

renormalizable operators characterized by mass scale A > Gp . Naively, the leading con- 
tributions to Df will be those suppressed by the fewest powers of [GpA"^)'^: namely, from 
dimension-six operators contributing to a[j^ that interfere with the SM amplitude a)^^.- There 
is only one such operator 68| : 



^dim6 = j^unn^dn iH^eD^H + h.c. , (13) 

where c is a complex coefficient. H is the Higgs doublet, D^^ is its covariant derivative, e 
is the antisymmetric tensor (ei2 = 1), and ^ denotes transpose acting on SU{2)l indices. 
Setting the Higgs field equal to its vacuum expectation value, Eq. ( IT3|) generates a coupling 
of the W boson to the right-handed charge current u^'^^dR, shown in Fig. [2l Integrating 
out the W boson, we obtain (recall we set Vud = 1) 

=^dim6 = - -^{uR'y^dReLl^h'eL + UR'y''diidLl^UL) + h.c. (14) 
The operator of Eq. ( IT3l) generates at order (G'i?A^)~^ contributions to both a]^^ and km: 

For all models that contribute to Dt via Eq. fll3p . EDMs are correlated with Dt in an 
otherwise model-independent way: 

\dn\ = 1 X 10"^^ ecm x (16a) 

\dHg\ = 7 X 10-^"^ ecm x\Dt/K\ (16b) 

Ic/dI = 4.5 X 10-2° ecm X lA/ft^l • (16c) 
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FIG. 2: (a) Effective UR'j'^djiW^ vertex arising beyond the SM: e.g., (h) left-right symmetric model 
with W-W' mixing; (c) exotic quarks ur, dn with non-standard SU{2)l x U{1)y gauge couplings 
that mix with SM quarks ur, dn; (d) R-parity violating MSSM with baryon number violation and 
squark left-right mixing. In each case, mixing insertions (involving the Higgs vev v) are denoted 
by (g). 



The current bound \dn\ < 2.9 x lO^^^ecm 22[ implies |-Dj/k| < 3 x 10^^, far below present 
sensitivities. 

This indirect limit on Dt applies to the following models: 

• Left-right symmetric models with a W boson that mixes with the W and couples to 
the right-handed quark charge current. 

• Models with exotic fermions with non-standard gauge quantum numbers, e.g., exotic 
S't/(2)i-doublet vector quarks u and d that mix with the usual u and d quarks. 



• -R-parity violating (RPV) MSSM with baryon number violation, described below |69| . 

The relevant diagrams are shown in Fig. |2l The first two models were studied previously in 
connection with D in Refs. [gI, Issl 34 1, and we do not describe them here. 



The RPV MSSM is defined by adding to the MSSM superpotential gauge-invariant and 
renormalizable terms that violate either baryon or lepton number (but not both, to avoid 
proton decay) j30|. Contributions to Dt are generated by the baryon number- violating 

WnPv = K,kU^Dpi, (17) 

where U^, are superfields corresponding to the (charge-conjugate) and o?-^ quarks 
of generation i,j, respectively. Shown in Fig. HI the leading contributions to Df arise at 
one-loop from diagrams involving third generation squarks t^^ji and bi^R- This contribution 
relies on mixing between gauge eigenstates, described by (see, e.g., Ref. fn\ ) 

.S:fmix = —'nT't {At sin f3 -\- fi cos (3) Vj^tR — mf, {Al cos f3 -\- fi sin (3) b^^^ h.c. (18) 

where tan /3 is the ratio between up- and down-type Higgs vacuum expectation values, Atj, 
and /i are MSSM mass parameters, and rrit {mi,) is the top (bottom) quark mass. For 
tan /3 ^ 1, we have 

V _ yus^suVtb mmb tan f3 n At 

- 2A^~ ■ ^^^^ 



Lepton number- violating terms have been studied previously in eonncction with the R coefhcient 70 1 
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assuming degenerate squarks with mass rriq and treating Eq. ( !T8|) perturbatively by mass 
insertion. Bounds on n-n oscillations constrain IA312I ^ 10~^ if all squarks have mass 



rriq = 200 GeV |72|. but this bound is relaxed if only third generation squarks are light; IA123I 



is unconstrained [30|. In principle Eq. (IT^ can be as large as C(10~^) for mq,At,fi ~ 200 
GeV, A" ~ 1, and tan/3 ~ 50 (perturbativity of the bottom Yukawa coupling requires 
tan/3 < 60). However, the neutron EDM bound constrains Im(a^^) < 3 x 10"'^, as per our 
previous discussion. 

Ref. [tsI previously studied the RPV MSSM in connection with EDMs, focusing on 
contributions from quark and electron CP-odd dipole moments arising at two-loop. For 
the combination of RPV couplings A" in Eq. f|T9|) entering Dt, quark EDM and chromo- 
EDM operators are suppressed by rriu^d- Here, the CP-odd four-quark operator gives a much 
stronger bound. 



IV. LEPTOQUARK MODELS 

Leptoquarks (LQs), fractionally-charged colored states carrying baryon and lepton num- 



ber, arise in many extensions of the SM, e.g., grand unification [7^ and compositeness [75|. 



Here, we consider a phenomenological model of LQs coupled to first generation quarks and 



leptons [32[. LQ models have a rich phenomenology for /3-decay, potentially giving large 
contributions to D and other observables through tree-level processes [6|. 

In the literature, LQ models have been regarded as an "EDM safe" source of CP violation 
that might generate D as large as present experimental limits, without fine tuning j6| . These 
previously considered models (dubbed the "usual scenarios") rely on LQ mixing to generate 
a dimension-eight operator contributing to a]^^ at tree-level, which interferes with the SM 
amplitude a[]^ jl]. In addition, scenarios involving LQs coupled to right-handed neutrinos 
can also generate via the interference of two new physics amplitudes a)^^. ^Yir- 

In this section, we study in detail these cases (i.e., with or without right-handed neutri- 
nos). We show that radiative corrections generate contributions to EDMs (in the spirit of 



Refs. [76|, l77|) sensitive to the same phases entering Dt. In both cases, the resulting bounds 



from the neutron EDM are stronger than the direct experimental limit. 



A. Usual LQ scenarios: no right-handed neutrinos 

There are two cases to consider: scalar and vector LQ exchange, both considered previ- 
ously in Ref. [cj. Since both cases are similar, we treat them simultaneously. The relevant 
LQs are 

scalar case: R= ( ) ~ (3, 2, 7/6) R= ( ) ~ (3, 2, 1/6) (20a) 
vector case: V = (^^^ ~ (3, 2, 5/6) \/ = ^+ ^ ~ (3, 2, -1/6) (20b) 
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FIG. 3: Scalar LQ case: tree-level exchange generates P-decay amplitude (left), while Olr is 
generated by one-loop vertex corrections (right), contributing to EDMs dn, dng, do- Diagrams are 
shown in weak-eigenstate LQ basis to illustrate that the same CP-violating phases from LQ-mixing 
( denoted fE> ) and couplings enter both Dt and EDMs. 

where ± states are weak isospin components, and the SU(3)c x SU{2)l x U{1)y quan- 
tum numbers are given in parentheses.^ In both cases, the most general renormahzable 
interactions to first generation fermions (including z/e/j) are 

scalar: ifjnt = URlleR + hp. QlCrR + hp dnLleR + Iir QLi^eRR + h-c (21a) 

vector: ^int = d^YLleV/, + gR Ql^f^eR eV^ + gi u'pfi^LleV^ + gR Ql^f^VeR eV^ + h.c. 

(21b) 

with quark and lepton doublets Ql = {uL,dL) and Ll = {I'eL^^Lj- Here, /iL,ij, gL,R-i 

gi^R are couplings (with L,R denoting lepton chirality). The presence of both L, i?-type 
couplings will lead to lepton universality violation in vr"*" e^v\ to avoid this constraint, we 
set i?-type couplings to zero j32|. The relevant mass terms are 

scalar: -=Sfmass = m\R^ R + m\R^R + {\r{R^ H) (RH) + h.c.) (22a) 

vector: ^^^ss = m^yViV^ + mlV^V^ + [\v{VIH){V^H) + h.c.) (22b) 

Through electroweak symmetry breaking, the quartic interactions (with couplings \Ry) give 
rise to i?_--R+ mixing and V^-V^ mixing by generating off-diagonal mass terms proportional 
to XRyv"^, where v = (H^). Diagonalizing the R--R+ and V--V+ mass matrices, we can 
express the mass eigenstates, denoted T^i 2 and Vi^2, as 

scalar: 7^l = cos Or i?„ + sin Or e''^«^+ , 7^2 = cos Or R+ - sin Or e~^*«i?_ (23a) 
vector: Vi = cos Oy V_ + sin Oy e'"^^ V+ , V2 = cos 9y V+ - sin Oy e-'^V_ (23b) 



^ We follow the notation of Ref. [32| for LQ states, except we omit an additional subscript identifying the 
SU{2)l representation. 
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FIG. 4: Vector LQ case: Tree-level exchange generates (3-decay amplitude (left), while Olr is 
generated by one-loop vertex corrections (right), contributing to EDMs dn, dug, do- Diagrams are 
shown in weak-eigenstate LQ basis to illustrate that the same CP-violating phases from LQ-mixing 
( denoted ® ) and couplings enter both Dt and EDMs. 

with mixing angles Buy and mass eigenvalues given by 

2\Xr\v^ 



scalar: tan29R = rnl^,, = l (ml + m| ± J {ml - + 4\XR\hA (24a) 

« R 

vector: tan 2dv = ""^Vi 2 = ^ {j^v + ''^1 ^ \J ("^v ~ + 4:\\v\'^v^ ^ 



2\Xv\v^ 

"^Vi,2 = 2 [my + m~ ± ^ [m^y - mt}' + AlXvl'V" j 

(24b) 

and phases = arg(ARy), defined such that m^y^ < rra^y^. The remaining (unmixed) 
LQ states R+, V+ and R-, V- have masses rriRy and "^5y, respectively. 
For /9-decay, this model gives i^t = Klm(aYji), where p] 

^ hLhlsm29ne^^^ (I I \ 
scalar case: Orr, = ^ — 7, ^ (25a) 

The relevant Feynman diagrams are shown in Figs. [3] and HI 

Next, we consider implications for EDMs. Radiative corrections involving the W boson, 
shown in Figs. [3] and HI generate the CP-odd four-quark operator Olr given in Eq. ([7]) 
which contributes to dn and rfng. The resulting coefficient kiR is proportional to the same 
CP-violating phases in Eqs. (!25l) entering D. For each case, we find 

SGpTTll / V \ SGpTTl^ / V \ 

scalar: kiR = —p^ FRlmfarp) , vector: kiR = 1= Fylmfarp) . (26) 

V2(47r)2 " ^ ^ V2{47iy ^ ^ ^ ' 

The loop functions F^y are given by 

= or^2 ^'2 X {fi"^n^ . rnl^ . "^r) + fi'^^n^ > "^Ti^ > "^|) + /("^Tii , "^7^2 . 0)) (27a) 
= Of 2 2 \ (3/("^Vi'"^V2'"^v) +3/(my^,mv2,m~) - /(m^^, m^^, 0)) (27b) 
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where 



£( 2 2 2\ 



log(m^/m|) + m^m^ Xo^im^jm 



m|)(m2 



mo 



(2^ 



Defined in tliis way, we have F^y > 1, with equahty in the hmit mj^y = 'ti-- ^ )^r,vv^ 



Eq. ( 1261) provides the leading contributions to EDMs from CP violation entering D; there 
is no suppression by light quark masses ^ as previously argued ^ . Using the results of 
Sec. [ni we have (both cases give the same numerical factors) 



\dJ > 4 X 10-21 ecm x \Dt/K\ ( """f Y 
' ' ' ^ ' VSOOGeV/ 

I^Hgl > 3 X 10"^^ ecm X \Dt/ k\ 
Idol > 1.7 X IQ-^i ecm X \Dt/n\ 



300 GeV 



300 GeV 



(29a) 
(29b) 
(29c) 



where miQ = rn-ji^ ("^vj corresponds to the lightest LQ state entering /3-decay for the scalar 
(vector) LQ case, (k ~ 0.87, —1.03 for rj^^Ne, respectively.) 

Recent searches at hadron colliders [tsI, 79, 80, [8l| provide constraints on the mass of 
the lightest LQ (T^i, Vi) involved in /3-decay. These bounds depend on the branching ratio 
/3e = BR(LQ ^ je) = 1 - BR(LQ ^ jv) = sin^ Ony, where j is a jet. For the scalar case, 
the strongest limits have been obtained at the Large Hadron Collider by combining jjee 
and jjeu channels (tsI . 79|: 



scalar case: m-n-^ > 



340 GeV 
319 GeV 



(CMS) 
(ATLAS) 



{(3e > 0.5) 



(30) 



with stronger limits (384 and 376 GeV, respectively) for /3e — )■ 1. Additionally, recent ATLAS 
searches for jets with missing energy from squark pair production, within a simplified SUSY 
context 12], apply to jjuu final states from TZi pair production. To translate the SUSY 
model into our framework, one must rescale the SUSY cross section by a factor (1 — /3e)2/4 
and take the gluino to be massive.'^ The resulting limits are m?^^ > 500 GeV, for /3e < 0.5, 
with stronger bounds in the limit /3e — ?■ 0. In the vector case, the DO collaboration has 
obtained 81, 83 



vector case: niy^ > 



302 GeV {jjee + jjeu, 0.1] 
144 GeV {jjuu, (3, < 0.1) 



(31) 



with stronger bounds for /3e — 1 or for different choices of anomalous gluon-LQ couplings 
considered therein. Within the context of our model, for (3e = sin^ < 0.1, the lightest 
vector LQ Vi ~ V- is nearly degenerate with 1+ (mv^ ~ my). Since BR(V"+ — > je) = 1, we 
have my > 367 8l|, and therefore is constrained indirectly to be much heavier than 
144 GeV. Additional constraints have been obtained by the HI collaboration at HERA. ^ 



^ The factor 4 counts the number of first and second generation squarks in the simplified SUSY model 
considered in Ref. [82[. 

we note S^/^ = R, S^j^ = R, V^i^ = V, 



^ To translate between the notation used here and that in Ref. 
and V^^i^ = V. 



11 




FIG. 5: Dt = K Im{a^^a^*^) is generated by LQ couplings involving right-handed neutrinos, with 
"iJL from R,R- or V,V -exchange (mixing denoted (^), and aj^^ from S- or U -exchange. 



These limits depend on the LQ-fermion couphngs, and provide stronger bounds than those 



from hadron coUiders if the relevant e-g-LQ coupling is larger than ~ few x 10~^ 84 . 

Given the current limit \dn\ < 2.9 x 10~^® ecm |22|, and conservatively taking m^Q > 300 
GeV, we conclude that \Dt/K\ < 7 x 10~^. CP violation in LQ models cannot saturate 
present experimental sensitivities in Dn — unless the hadronic uncertainties associated with 



the dn computation of Ref. [35[ are larger by an order of magnitude, or unless there is 
a cancellation with other CP-odd phases contributing to dn to ~ 10% (or a combination 
thereof). On the other hand, the mercury EDM does not strongly constrain Dt in this 
model, especially in light of its large hadronic uncertainties, although this situation may 
change with future refinements in the nuclear computations. 



B. LQ scenarios with right-handed neutrinos 

LQ models can contribute to Dt through the interference between two new physics am- 
plitudes involving right-handed neutrinos. The relevant Feynman diagrams are shown in 
Fig. [5l To begin, we consider the model of the preceding section involving scalars i?, R and 
vectors V, with mixing defined in Eqs. fl22ti24p and couplings to SM fermions given in 
Eq. (!2T|) . Here, we set to zero L-type couplings in Eq. ( 12T|) and keep only i?-type ones. For 
each case, the amplitude a)^^ is 

^ , V sin 2e^e-^'^« / I 1 \ , , 

Tci 2-exchange: a^r = — — ^ ^ (32a) 

V.,-exchange: a,, = ^^^^ " " ^'''^ 

In order to generate ct/j^, we introduce two additional LQ states S and f/, with quantum 
numbers 

scalar LQ: S ~ (3, 1, 1/3) , vector LQ: U ~ (3, 1, 2/3) (33) 
and quark-lepton couplings 

=^int = {gs UrCr + g's d^^Ueu) S + [hu dnl^CR + h'^j unY^eR) + h.c. (34) 

Through tree-level exchange, these states generate 

S-exchange: = -^i-^ . [/-exchange: = -_^L^ . (35) 
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FIG. 6: LQ contributions to Dt generate radiatively CP- odd operator Olr contributing to EDMs. 



There are four possible contributions to Dt = fi:Im(a^^a)^|j) depending on which of the 
combinations 

(7^l,7^2,^), (7^l,7^2,f/), (Vi,V2,^), {v,,V2,u) . (36) 

we consider contributing to a^j^ and o^r- 

Next, we consider each of these combinations separately and compute the resulting EDM 
induced by the CP-odd four quark operator in Eq. ([7]). There are four possible contributions, 
shown in Fig. [6l and they all give nearly identical results: 



I y^GprnlQ I , V v*\\ ft 2 2 ■ 



(47r) 



(37) 



The loop function is 



f{ml,ml,ml 



2 rn\rn2rn\{m\ log(m2/m3) + m\ \og{m\/mi) + m\ log(m^/m2)) 
Tn\Q{m\ — mf)(m3 — m^{m\ — m|) 



(38) 



where, for each case in Eq. f p6|) . ml 2 3 corresponds to the masses of the three states, with 

m^Q being the smallest of the three. Defined in the way, we have / > 1, with equality if all 
states are degenerate. 

Assuming that one CP- violating phase is dominant in Dt, the resulting EDMs arising 
from that phase are 



M„| > 9 X 10- ecn X |A/K| (^^)^ 
|d„,|>7xlO-ecmx|A/K| (5^)' 
|.„|>4xlO-ec,„x|A/«| (3^)' 



(39a) 
(39b) 
(39c) 



Comparing Eqs. fl26l) and (137|) . we find that Dt from LQ scenarios involving right-handed 
neutrinos is less constrained by EDMs by a factor 4 compared those involving left-handed 
neutrinos (for fixed TTi^g). 

Constraints on scalar and vector LQ masses from pair production at hadron colliders are 
the same as in Eqs. ( l30l) and ( 13T]) . However, in the limit hu <^ h'^, the vector U decays 



primarily via U — )■ ju and is subject to the relatively weaker mass bound rriu > 144 GeV [83 
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Significantly stronger bounds are provided by the HI collaboration for (3e{U) ~ 8J|, which 
depend on the U-e-d coupling hu- 



( 250 GeV {hu = 0.03) 
mu > < 300 GeV {hu = 0.06) (40) 
[ 1 TeV {hu = 0.3) 

Although suppressing hu weakens the bound on mu, the contribution to Dt{(x hu/mfj) 
is also suppressed. Assuming hu > C(0.06) (to avoid too much additional suppression 
in Dt) we take mu > 300 GeV.^ For miq > 300 GeV, the neutron EDM bound implies 
|A/«:| <3 X 10-5. 



V. CONCLUSIONS 



The emiT collaboration has measured = (—1.0 ± 2.1) x 10 [15j, consistent the 
SM prediction dominated by O{10~^) final state effects. Here, we studied several new 
physics scenarios beyond the SM and showed that the current neutron EDM measurement 
\dn\ < 2.9 X lO'^^ecm provided in all cases stronger bounds on D. 

• l-Df/^l < 3 X 10"'' in left-right symmetric models, exotic fermion models, and the 
R-parity violating MSSM. EDM bounds on this class of models, given in Eq. f|T6l) . 
can be understood in an otherwise model-independent operator framework through a 
coupling of the W boson to the right-handed quark charge current UR'^^^dR. 

• \Dt/ n\ < 3 X 10"^ (7 X 10^^) in leptoquark models with (without) light right-handed 
neutrinos. Moreover, EDM constraints will become more severe if collider bounds on 
leptoquark masses are improved, as shown in Eqs. fl29|) and (139]) . 



We recall that k ~ 0.87 (for the neutron) is defined in Eq. ([2]), and Dt denotes the contri- 
bution to D from fundamental T violation (as opposed to final state effects). Analogous 
constraints from the mercury EDM bound are weaker by an order of magnitude (with large 
uncertainties), although the situation may change with future improvements in the nuclear 
computations. A future constraint on the deuteron EDM of \dD\ ^ lO'^^ecm would im- 
prove all aforementioned bounds on Dt by two orders of magnitude. These bounds can in 
principle be evaded by fine-tuned cancellations with other CP-odd phases contributing to 
EDMs, but not to A- 

Even though D is not as sensitive as EDMs to CP violation beyond the SM, clearly it 
worthwhile to push D measurements to greater sensitivities. Since any single EDM measure- 
ment has little model discriminating power, it is desirable to consider as many observables 
as possible — especially if a non-zero EDM were measured. In this case, D could play an 
important role in untangling the nature of CP violation and potentially shedding light on 
origin of matter in the Universe. 



^ It seems plausible that the best trade-off between small mu and small hu oecurs for rau ^ 300 GeV, cor- 
responding to the center-of-mass energy = 319 GeV at HERA. For rajj < y/s. on-shell LQ production 
dominates, allowing for relatively stron ger constraints on hu] for mjj > ^/s, only off-shell production is 



allowed, and the constraints arc weaker 8^. A more precise analysis is beyond the scope of this work. 
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